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Abstract. A degenerate nonlinear nonlocal evolution equation is considered; it can be understood as a porous 
' medium equation whose pressure law is nonlinear and nonlocal. We show the existence of sign changing weak 

i solutions to the corresponding Cauchy problem. Moreover, we construct explicit compactly supported self- 

similar solutions which generalize Barcnblatt profiles — the well-known solutions of the classical porous medium 
1 equation. 
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1. Introduction 

| In this work, we study the following degenerate nonlinear nonlocal evolution equation 



> 

On 



(1.1) d t u = \7 ■ (\u\V a - 1 (\u\ rn - 2 u)) , xeR d ,t>0, 

where m > 1 and V a_1 denotes the intcgro-diffcrcntial operator V(— A)^ _1 , a £ (0,2). The equation is 
supplemented with an initial condition 

(1.2) u(0,x) =u (x). 

£ — , First, we construct nonncgativc self-similar solutions of equation (jl.ll) which are explicit and compactly sup- 

ported. They generalize the classical Barenblatt-Kompaneets-Pattle-Zel'dovich solutions of the porous medium 
equation, see (|1.4p below. Second, we prove the existence of sign changing weak solutions to problem (|l.ip - (|1.2p 
for merely integrable initial data, and we prove that these solutions satisfy sharp hypercontractivity L 1 i— > LP 
, estimates. 



A nonlocal operator. Equation (|1.1[) involves a nonlocal operator denoted by V a_1 which can be defined 
as the Fourier multiplier whose symbol is «£|£| a ~ 2 - This notation emphasizes that it is a (pseudo-differential) 
operator of order a — 1. Recalling the definition of the fractional Laplace operator (— A)% (v) = .F _1 (|£| a .F?;) 
and the Riesz potential Ip = (— A) - '? , i.e. Fourier multipliers whose symbols are |£| a and |£| -/3 respectively (see 
for instance [231 Ch. V]), the fractional gradient V Q_1 can also be written as V22-a- Finally, let us emphasize 
that the definition of V Q_1 is consistent with the usual gradient: V 1 = V; the components of V° are the Riesz 
transforms; moreover we have V ■ V" -1 = V' ■ = —(—A)'. It is also possible, following the reasoning from 
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[12\ Th. 1] to define the fractional gradient via the singular integral formula for smooth and bounded functions 
v :R d ->■ R 

(1.3) V«-^(x) = C d , a J (v(x) - v(x + z)) dz 
with a suitable constant Cd, a > 0. 

Related equations and results. Our preliminary results for the problem (|1.1[) - (|1.2[) have been announced 
in 0. 

First, we would like to shed light on the link between and other partial differential equations. Notice 
that when a = 2 equation coincides with the classical (nonlinear parabolic) porous medium equation 

(1.4) dtu = V- (\u\V(\u\ m - 2 u)) = V- ({m-l)\u\ m - 1 Vu). 

For the theory of porous media equations, the interested reader is referred to [37J HB] and references therein. Of 
course, for m = 2, the Boussinesq equation is recovered. 
The following nonlinear and nonlocal equation 

(1.5) 9iW + |t,l| (^) 2 " = ° 

in the one-dimensional case x G R was studied by the first, the third authors and R. Monneau [1]. Such an 
equation was derived as a model for the dynamics of dislocations in crystals. In [4], the existence, uniqueness 
and comparison properties of (viscosity) solutions have been proved, and explicit self-similar solutions have been 
constructed. Notice that the function u = v x , where v is a solution to (|1.5[) . solves the one-dimensional case of 
p. II) with m = 2. Thus, equation (jl.lj) is a multidimensional generalization of the one in (|1.5[) . 

Recently, Caffarelli and Vazquez [SJ [7J studied nonnegative weak solutions of (|1.1[) in the case m = 2 in the 
multidimensional case. Precisely, they studied the following (nonlocal) porous medium equation in R d 

(1.6) d t u = V ■ (uVp), 

with the nonlocal pressure law p = (— A)~ s u, < s < 1, obtained from the density u > 0. Notice that, 
for a = 2 — 2s G (0,2), equation (|1.6[) reads dtu = V • (mV" _1 u). For sign changing u's, our equation (|1.1|) 
is a (formally parabolic) extension of equation (|1.6[) of the structure of (|1.5[) . In [5], Caffarelli and Vazquez 
constructed nonnegative weak solutions for (|1.6[) . i.e. for (jl.l[) with m = 2, with initial data satisfying: ito G 
i 1 (R d ) fl L°°(R d ) and such that < Uo(x) < Ac~ a ^ x ^ for some A, a > 0. Besides the positivity and the mass 
preservation, the properties of solutions, listed in the next paper [7J p. 4], include the finite speed of propagation 
proved using the comparison with suitable supersolutions. Further regularity properties of solutions of Ql.ip 
with ui == 2 and a G (0, 1) are studied in [B]. 

Another nonlocal porous medium equation has been proposed in |10l II H 129] 

d t u+ (-A)t(| M | m - 1 u ) = 



NONLOCAL POROUS MEDIUM EQUATION 3 

for a = 1 and a G (0, 2), respectively. Among several other properties like smoothing effects and decay estimates, 
solutions of this generalization of the porous medium equation enjoy the L -contraction property, so, they are 
unique. But self-similar solutions are not compactly supported [291 Th. 1.1]. 

Finally, we recall that the following nonlocal higher order equation, appearing in the modeling of propagation 
of fractures in rocks, 

d t u = V- (u' l V(-A)5 u ) 

(with u > and n > 1), has been studied in |15j in a one-dimensional bounded domain. At least formally, this 
equation with n = 1 corresponds to (|1.1[) with a = 3 and m = 2. 

Notation. In this work, Q T denotes (0,T) x R d . The usual norm of the Lebesgue space L p (M. d ) is denoted 
by || • Hp for any p G [l,oo], and H s ' p (M. d ) with the norm || • \\h°-p is the fractional order Sobolev space, see 
Section [3] The Fourier transform J- and its inverse transform T~ x of a function v G are defined by 

Tv(0 = (2tt)-* J v(x)e- ix£ dx, F- X v{x) = (2ir)~i J v(x)e tx< d£. 

Here, all integrals with no integration limits are over the whole space R d if one integrates with respect to x 
and over the whole half-line R + = [0, oo) if the integration is with respect to t. As usual, w+ = max{0, w}, 
W- = max{0, — w}, sow = w+ — W-. Constants (always independent of x and t) will be denoted by the same 
letter C, even if they may vary from line to line. Sometimes we write, e.g., C = C(p,q,r) when we want to 
emphasize the dependence of C on particular parameters p, q, r, for instance. 

2. Main results 

In this work, we show two main results: we construct explicit self-similar solutions of equation (jl.ll) . as well 
as we prove that the initial value problem (|l.ll) - (jl.2l) has a global-in-time weak solution which satisfies certain 
optimal decay estimates. 

We first "recall" the appropriate notion of weak solutions for Equation (jl.ip , see for instance [27j [28] . 

Definition 2.1 (Weak solutions). A function u : Qt — ^ K is a weak solution of the problem (|1.1[) - (|1.2[) in Qt 

Hue L^Qt), V Q - 1 (|u| m - 2 u) G L\ oc {Q T ) and \u\V a - 1 (\u\ m - 2 u) G L£, c (Qt). and 

J J (udttp- \u\V a - 1 (\u\ m - 2 u) ■ Vtp) dtdx + J uo(x)<p(0,x)dx = 

for all test functions ip G C°°(Qt) H C{Qt) such that <p has a compact support in the space variable x and 
vanishes near t = T . 

The first main result of this work says that there is a family of nonnegative explicit compactly supported 
self-similar solutions of (|1.1|) . i.e. nonnegative solutions that are invariant under a suitable scaling. Observe 
that if u(t, x) is a solution of ([l.ljl . then so is L dx u(Lt, L x x) for each L > 0, where A = (d(m — 1) + a) 1 - Thus, 
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the scale invariant solutions should be of the following form 
(2.1) u ( t ,x) = -^$(£) with A 



t dx \t x J d(m-l) + a' 

for some function $ : R d — > M. satisfying the following nonlocal "elliptic type" equation 

(2.2) -AV-(i/$) = V- (I^IV"-^!*!'"- 2 *)) where y = ^. 

Theorem 2.2 (Self-similar solutions). Let a E (0,2], m > 1. Consider the function & a ,m : K d — > R defined as 

(2.3) $a,m(lO = (kad 1 ~ \V\ 2 ) + 



with the constant 

dT(l) 



(d{m - 1) + a)2 a T (l + f ) T (*t2) ' 
Then, the function u : (0, oo) x W l — > R + defined by (|2.1[) $ = $ Qi , n is a weafc solution of (jl.lj) m i/ie sense 

of Definition\MJ\ in Q V . T = (r), T) x K d /or every < n < T < oo. Moreover, u(t,x) satisfies the equation in the 

pointwise sense for \x\ ^ t d ( m - 1 )+» , and is min { 2(w"-i) ' lf-Holder continuous at the interface \x\ = i)+a . 



Remark 2.3. When a = 2 in expression (|2.4[) below, we recover the classical Barenblatt-Kompaneets-Pattle- 
Zcl'dovich solutions of the porous medium equation (|1.4[) , see for instance [28l [27] . 



Remark 2.4. For each M E (0, cx>) we can find a nonnegative self-similar solution u with prescribed mass 
M = J u(t, x) dx (which is conserved in time) by a suitable scaling of the profile & a ,m- Indeed, this self-similar 
solution is given by the formula 



(2.4) u(t,x) = t -Hm-D+c. ( k a d ( e 2 - Xt Hm-lU 

where, for each M > 0, there exists a unique R > such that J u(t, x) dx = M. 

Remark 2.5. Self-similar solutions of equation (|1.6|) (which is a particular case of equation Ijl.ljl ) have been 
proved to exist in [7] by studying the following obstacle problem for the fractional Laplacian. For a E (0, 2) and 
^(y) = C — a\y\ 2 where a = a(d, a) and C > 0, one looks for a function P = P(y) with the following properties: 

P>*, (-A)tp>0, and either P = f or (-A)*P = 0. 

The novelty of our approach is that we exhibit the explicit self-similar profile $ a 2 defined in (|2.3[) and, conse- 
quently, the explicit solution of this obstacle problem: P(y) ~ T a ($0.2) (-f|) , where I a = (— A)~t is the Riesz 
potential and R > is a suitable constant. 

Next, we prove the existence of weak solutions to the initial value problem (|1.1[) - ([1.2[) . 

Theorem 2.6 (Existence and decay of L p -norms). Let a E (0,2) and 

' m> 1 + 1=2 i/ a£ (0,l], 



(2.5) 



m>3-§ ifae(l,2). 
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Given uq £ L 1 (M d ), there exists a global-in-time weak solution u of the Cauchy problem (|l.lj) - (|1.2[) . Moreover, 
and 

(2.6) IK*)IIp < C(d,a,77i)||uo|li^ m_)l ''° r '<™-V» C 1 "?) for all t > 0, 

ZioZds with the constant C{d,a,m) independent of p and uq- 

The solution u is nonnegative if the initial condition uo is so. If uq £ L p (M. d ) for some p £ [1, oo], then 

K*)IIp<II«oIIp 

holds for all t > 0. 

Remark 2.7. Estimates (|2.6p are sharp since the decay in Theoreni l2.6l corresponds exactly to that for self-similar 
solutions constructed in Theorem 12.21 Moreover, for a = 2, they are similar to those for degenerate partial 
differential equations like the porous medium equation (showing the regularization effect on the L p -norms of 
solutions); see, e.g., [23], [SJ Ch. 2]. 

Remark 2.8. After proving those hypercontractivity estimates in [3J, we learned that a similar result is obtained 
in [BJ, however, for a less general model: a £ (0,1), m — 2, and nonnegative uq. Moreover, analogous decay 
estimates for another fractional porous medium equation of the form dtu + (—A)%(\u\ m ~ 1 u) = were proved 
recently in 

Compared with the methods used in [3] , we propose an alternative strategy of the proof of the existence of 
solutions. In this paper, we consider approximating solutions u = u S ' e of the equation 

dtu = 5Au + V- (|it|V Q_1 G e (u)), 

considered in the whole space R d , where G e (u) is a sufficiently smooth approximation of u\u\ m ~ 2 , and then we 
pass to the limit with the parameters e \ 0, and S \ 0. Solutions of the approximating equation exist because 
the parabolic regularization term SAu is strong enough to regularize equation when < a < 2, but of 
course not for a — 2. Our approach resembles the approach to the one-dimensional model achieved in [?1 Sec. 
4 and 5] via viscosity solutions. 

3. Preliminaries 

In this section, we collect known results that we will used in proofs of the main theorems. 

Bessel and hypergeometric functions. Bess el functions of order v are denoted by J v (z), and they behave 
for small and large values of the (complex) variable z like 

J v {z) - TW+T)(l) U as z "> °. 

Jv(z) ~ (f)5 cos(z - t£ - j )z~5 as \z\ -)• oo, 
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where for functions /, g, the relation / ~ g means that ^ — > 1. For the proofs of those properties of J„, the 
reader is referred to, e.g., |30j . 

The hypergeometric function, denoted by 2^1(0, 6; c; z), is defined for complex numbers a,b,c and z as the 
sum of the series 

2 F 1 (a,b;c;z) = f^^^z" for \z\ < 1, 

where (a)„ = ^^y^, and T denotes the Euler Gamma function. This series is absolutely convergent in the 
open unit disc and also on the circle \z\ — 1 if 5i(a + b — c) < 0. 

It is known [19[ p. 39] that when b = —n is a negative integer, 2-F1 (a, —n;c;z) is a polynomial function of 
degree n. In particular, we have 

(3.1) 2 F 1 (a,-l;c;z) = l--z. 



c 



We will also use the following differentiation formula [TH1 p. 41] 



(3.2) ^(aFi(a,&;c;*)) = ^ 2 Fi(a+ 1,6 + l;c+ l;z). 

77ie Weber-Schafheitlin integral. If < 6 < a and if integral (|3.3p below is convergent, then the following 
identity holds true 

f°° x T/ w/ , b^2- A a A - t/ - 1 r( t/+ ^7 A+1 ) /l + u-A + I^-u-A + I 6 2 

(3.3) / t- x J ll (at)Ubt)dt= ^ r _„ +u+x+1 ^ 2 , V i — ^- , — ^- ;^ + i;-n 



r( -»+H-*+ 1 )r(i + i/) V 2 ' 2 

According to Watson, [301 PP- 401-403], this result was obtained by Sonine and Schafheitlin. However, it is 
usually referred to as the Weber-Schafheitlin discontinuous integral since there occurs a discontinuity for a = b. 

The Stroock Varopoulos inequality. We next recall the the Stroock-Varopoulos inequality, see [HI Theo- 
rem 2.1 and Condition (1.7)] for a proof. 

Proposition 3.1. For a 6 (0,2], w g C£°(IR d ) and g > 1, the following inequality holds true 



(3.4) J sgnw M 9_1 (-A)* wdx > 4(g ^ 2 ^ ^ |v* (sgnu> |«;|*) |* dx > / 



V 2 M 2 



dx. 



Fractional order Sobolev spaces. The fractional order Sobolev spaces are defined as 

H s < p (R d ) = {v E L p {R d ) : V s v G L p (R d )} = {»£ L p (R d ) : (I - A)^v € L p (R d )}, 

here with p £ (1, 00), supplemented with the usual norm denoted by || • ||jjs,p, and we refer the reader to the 
books [23 EB] for properties of those spaces. In particular for s = a — 1 with a G (1, 2), the following well-known 
continuous embedding will be used repeatedly 

(3.5) H a - Up (R d ) C L°°(R d ) provided p > — ^— (> 1). 

a — 1 

We also recall the fractional integration theorem [23j Ch. V, §1.2]: the Riesz potential I s = (— A)~i satisfies 

(3.6) \\2 s u\\ q <C(p,q,s)\\u\\ p 
for all s e (0, d) and p, q € (1, 00) satisfying | = i — ^ . 
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Some functional inequalities. We will use the following Nash inequality 

(3.7) Hl2 (l+f) < cwiiv^miH? 1 

valid for all functions v £ L 1 (R d ), such that V^w £ L 2 (R d ), and with a constant Cn = C(d, a) > 0. The proof 
of (|3.7|) for d = 1 can be found in, e.g., jTHl Lemma 2.2], and this extends easily to the general case d > 1. 
Moreover, we will use the following Gagliardo-Nirenberg type inequality 

Lemma 3.2. For p > 1 and p >m—l, the inequality 

(3.8) IK<^||V*H5||"||u||5 
holds with 

p d(r-l) + a d(m- 1) + pa 

(3.9 a = , b = a-r= — , r = p + m-l. 

p — 1 d d{p — 1) 

Proof. This inequality is a consequence of the Nash inequality (|3.7p written for v = |it|3 , i.e. 

(3.10) <C N \\V^\u\m 2 2 \\u\\f, 

and two Holder inequalities 

M^wurjuw 1 ^ with 7 

and 

W*<ll«||jHll- a with 6 

* \p — 1 j r 

Combining the above three inequalities, we get (|3.8I) . □ 

4. Proof of Theorem 12.21 

This section is devoted to the study of nonncgative self-similar solutions for with m > 1. As explained 
above, this problem reduces to a study of the elliptic-like equation (|2.2j) which for nonncgative $ takes the form 

Moreover, since we want to construct compactly supported solutions, we are interested in solutions $ vanishing 
outside the unit ball B\ . This is the reason why we consider the Dirichlet problem 

-Ay = V"- 1 ^™" 1 ) in B x , 

(4.1) 

$ = in K d \ Si . 

It is well known that, in the case of nonlocal operators (such as V a_1 ), the homogeneous Dirichlet condition 
should be understood in the form $ = outside the domain B\, and not only $ = on the boundary dB±. 
The reader is referred to, e.g., [5] for more explanations. 

We claim that the proof of Theorem 12.21 reduces to the following key computation. Here, 2F1 denotes the 
classical hypergeometric function defined in Section [3] 



r — 1 J p' 
r — 2 \ p 



(4.2) l p (l-M 2 )! = 
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Lemma 4.1. For all (3 G (0, 2) , j3 < d, and 7 > 0, we have 

C 7Ad X 2 ^(^,-2M;|;|y|2) for | y | < ^ 

C- 7 , W |yr d X 2 ^(^,^;^;^) for |y| > 1, 

r(| )r(£±^+i) aM °T^> d ~ z r (f)r(±±i+i) ' 
The proof of this lemma is postponed to the end of this section. 

The following corollary is an immediate consequence of Lemma |4 . 1 1 with 2 — /3 = a = 7, of the property of 2 i*i 
formulated in (|3.1j) . and of the identity (—A)" = (— A)X 2 -a- It has an important probabilistic interpretation, 
and recently, related results and generalizations have been proved in |13j . 

Corollary 4.2 (Gctoor [H Th. 5.2]). For all a G (0,2], i/ie identity 

K a , d (-A)%(l-\y\ 2 )l =1 in Bj 
irae wi/j/i f/ie constant K a ^ = — — 7 — g d+et ? . 

Before proving Lemma [4~T| we first use it to derive Theorem 12.21 

Proof of Theorem \2.SX We check that u(t,x) = t~ dx Q a . m (t~ x x) is a weak solution of (jl.ll) in the sense of 
Definition 12.11 First, u G L 1 (Qt) if and only if $ Qim G L 1 (K d ), which is obviously true. For later use, it is 
convenient to introduce the function $ Q (w) = (l — |y| 2 )^, so that k a ^a = "J*™™ 1 - 

The fact that, for all n,T such that < rj < T, uV*- 1 ^™- 1 ) and V a " 1 (u ro_1 ) are locally integrable in 
(77, T) x R d follows from 

which we prove by computing Z 2 -a:(^a:)- In order to do so, we first assume that a > 2 — d, and we apply 
Lemma H. II with 7 = a G (0, 2) and /? = 2 — a, we use equation ()3.1[) . and we get 

(4.3) 2a-a(*aX») - j C a , 2 _ Q , d | y | 2 -( d +«) 2 F X - 1, f ; ^; ^) if |y| > 1. 

The right-hand side of equation (|4.3[) defines a locally integrable function because 

,d + a . a d + a 
a + *>-c=(— -1) + --— <0. 

We then deduce that 

V Q - 1 ($ Q )(zj) = VX 2 _ Q ($ Q )( 2 /) = --^j/ for yGfli. 

Note also that V Q_1 (<I , Q ,) can be computed outside B\ thanks to the differentiation formula (|3.2[) . 
We now remark that 

$ a , m (y)(V a - 1 ^- 1 )( 2 /) = -Ay$ a , m (y) for all y£< 
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which is in L 1 (]R <i ). Moreover, the following equalities hold true in the sense of distributions in Qt, 

d t u(t,x) = -Xt-^Wy ■ (y<S> a , m )(t~ x x), 

This allows us to conclude that u is indeed a weak solution of in (77, T) x R d for all < 77 < T < 00 if 

a > 2 - d. 

Assume now that < a < 2 — d, which means that d = 1 and a < 1. The critical case a = 1 can be obtained 
by passing to the limit as a \ 1; indeed, the constants C a ^-a,d(d + a — 2) and C a ^-a,d(d + a — 2) appearing 
in (|4.3p simplify thanks to the relation zT(z) — T(z + 1). If now a < 1, we can argue as above by analytic 
continuation. The proof is now complete. □ 



Now we turn to the proof of the main technical lemma. 

Proof of Lemma \JH\ We first assume that /3 € (0, d) and that 7 > max{0, d — 2/3 — 1}, and we then argue by 
the analytic continuation with a choice of parameters corresponding each time to 2-F1 defined and bounded for 
all \y\ < 1. 

The Fourier transform of $ 7 (y) = (1 — |y| 2 )^ is expressed in terms of Bessel functions, see, e.g., [22 Ch. 
IV, Sec. 3] 

J-(4> 7 )(£) = 2lT (1 + l) * 

Since Ip is the Fourier multiplier of symbol ^3(^7) is the (inverse) Fourier transform of the following 

radially symmetric function 

2*r (1 + 1) — i J<+t(|£|). 

We recall that by properties of Bessel functions collected in Section [3l we have J„(r) = O (r v ) as r — > and 
Jyir) = O (r as r — > 00. We see that the previous function is intcgrable since j3 < d and d < 7 + 2j3 + 1. 
Thanks to H Th. 3.3], we get 

\ l / £ — HP 2 2 

= 2*r (2 + 1) \ v \'-i J^t-^j^J^tlyDdt. 

We obtain (|4.2|) applying p.3[) with the following choice of parameters: 

• if \y\ < 1, we put A = 2 + /3, /j = ^jp, 1/ = | — 1, a=l and 6 = 

• if \y\ > 1, we put A = \ + /3, /it = f - 1, 1/ = ^jp, a = and 6=1. 

□ 
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5. A REGULARIZED PROBLEM 

In order to construct weak solutions of for general initial data, we first consider the following regularized 
problem 

(5.1) ftu = <yAu + V-(MV a_1 (G(u))), u(0,x) = uo(x), 

where G : K — > K satisfies 

{G differentiable and increasing, 
G(0) = G'(0) = 0, 
G' locally Lipschitz continuous. 

Remark that for m > 3 or m = 2, the function G(u) = \u\ m ~ 2 u satisfies (|5.2j) . For m G (1, 3), we consider the 

following approximation G = G £ of |u| m ^ 2 u 

G e {u) = sgnu {{u 2 +e 2 ) S ^ 1 - e'"" 1 ^ 

with e > 0. The following theorem holds true for a general function G satisfying (|5 . 2[) . 

Theorem 5.1 (Existence of solutions to the regularized problem). Let S > and assume that G is an arbitrary 
function satisfying (|5.2p . Moreover, assume 



L 1 (R d ) n L°°(R d ) i/ ae(0,l], 

L 1 (M d )n(n p>PQ if Q - 1 -f(R d )) ae(l,2), 



(5.3) u G 

with p a = — 3j > 1 . There exists a unique function u in the space 

(5.4) u G 
satisfying problem (|5 . lj) m i/ie usual weak sense 

(5.5) / / (u^tv? - |u|V a_1 (G(u)) • Vy> - <5Vu ■ Vv?) didx = 



C ([0, oo),L 1 (R d ) n L°°(R d )) if a £ (0, 1], 

n P > PQ C([0,c3o),i 1 (IR d )ni? Q - 1 'P(R d )) i/ a 6(1,2), 



for all<peC™(Q T ). 

Moreover, u(t, x) is nonnegative if the initial condition uo is so, and for all t > and q G [1, oo] we have 
(5.6) y u(t,x)da;= / uo(x)dx and \\u(t)\\ q < \\uo\\ q - 

Local-in-time existence of mild solutions. 



C([0,T),L 1 (R d )nL co (R d )) if 6(0,1], 
C ([0, T), L 1 (M rf ) n i2' a - 1 'P(R d )) i/ a £ (1, 2) 



Proposition 5.2. Lei p > p a = There exists T > depending only on uq, and a function u in the space 

(5.7) u G 
smc/i £/ia£ 

(5.8) u(t)=e 5tA u + [ Ve s ^ A • *(u(*)) da with = MV""^^), 

Jo 

m C[0,T],L 1 {M d ) n L°°(IR <i )) w/iere e tA denotes the heat semigroup. 
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Remark 5.3. Wc identify the heat semigroup c* A and its kernel (4nt)~% cxp ^— ^-J ■ Wc will use the following 
classical fact 

(5.9) || V e stA v\\ p < C{ P , r, (3, a)H(MM||t;|| P 

with 1 < r < p < oo, and j3 G [1, 2). 



Now we turn to the proof of Proposition 15.21 
Proof of Provosition \5.2\ We look for a solution u G C([0, T], X) as a fixed point of the map 

T:ii4 c 5tA u + / Vc 5 (*- S ) A ■ *(«(«)) ds, 
Jo 

where X is chosen as follows 

x= fLH^ni 00 ^) if a g (0,1], 

1i x (m d ) r]H a -^P(R d ) if a G (1,2). 

The associated norms are \\u\\i + \\u\\ Y with Y = L°°(R d ) and Y = ff Q| - 1 ' p (R rf ), respectively. We show that T 

has a fixed point by the Banach contraction principle as soon as T = T(||tto||x) > is sufficiently small. 

In both cases, it is enough to prove the following lemma. 

Lemma 5.4. For all T G (0, 1), the operator T maps C([0, T], X) into itself. Moreover, there exist C > and 
7 > such that for all u,v G B(0, i?) C C([0,T],X), 

(5.H) ||7» - 7»|| c([oM < C^R^Wu - HIc([o,t],x), 

where C\(R) is a constant which also depends on a,d,e,m,S (and on p if a G (1,2)). 

Indeed, once this lemma is proved, we first derive 
(5-12) ||7»|| c([0 ,t],x) < \\u \\x + RC 1 {R)T 

by choosing v = in (|5.11j) and using estimate (|5.9j) . Now it is enough to choose R = 2 1 1 uq \ \ x and T > such 
that Ci(i?)T 7 < \ in order to ensure that T maps _B(0, i?) into itself, and is a contraction. 

The case a G (0, 1]. In order to get estimate (|5.11j) . we first write 

(5.13) 7»(t) - 7»(f) = f Vc j ^ s ) A • (*(«) - («)) (s) ds, 

Jo 

and the difference of 'S's is represented as 

(5.14) #(u) - V(v) = (\u\ - |«|)V a_1 G(ti) + |w|V Q - 1 (G(u) - G(v)). 

Lemma 5.5. For every a G (0,1] and p G (1, oo) i/iere exists a constant C(p,a) > such that for all u G 
L°°(R d )r\LP(M. d ) the following inequality 

(5.15) ||V«- 1 (G r («)-G(«))|| g <C(p J a)[ sup G'(z) ) ]|u -«|| p 

\M<IMU + |M|oo / 

Zio/cfe £rue wi£/i - = ± — i^ 2 . 

IP d 
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Proof. Since, for a G (0, 1], we have V Q_1 = V°Zi- Q , where the components of V° are the Riesz transforms 
(which are bounded operators on L p (M. d ) for each p £ (1, oo)), we obtain (|5 . 1 5[) from estimate p.6[) as follows 

HV"- 1 ^) _ G ( v ))\\ q < ||Z!_ a (G(tt) - G(v))\\ q < C( P , a)\\G(u) - G(v)\\ p <C\\u- v\\ p . 

□ 

Now, we come back to the proof of (f5"TTj) with X = L 1 (M d ) n L°°(K. d ), First, for all u, v G B(0, R) C X and 
some q G (1, oo), 

»*(«) - *(«)||i < - liv"- 1 ^)!!, + iiv-^Gfu) - Gf(«))||, 

(5.16) < C\\u - t'W^G'iWuW^WuWp + ciMi^c'dHio, + ||w||oo)H« - v\\ P 

< C (R)\\u-v\\ x 

with C (R) = C(a, d, q, RG'(2R)) and ± + ^ = 1 and i = | - i=2. We used estimate ([5~T5j) twice to get the 
second line in (|5.16[) . and the inequality ||u|| r < ||u||x which is valid for all r G [1, oo] to obtain the last one. 
The estimate of the second norm in X is obtained similarly: for all u, v G B(0, R) C X and some q G (1, oo), 

(5.17) ||*(u) - *(«)||, < ||u- wIUIIV -^^)!!, + MUIIV^^G^) - G(v))\\ q < C (R)\\u - . 

Now, we apply inequalities (|5 .9[) with /3 = 1, (p, r) = (1,1) and (p, r) = {oo,q), respectively, and the estimates 
(15161) . (15371) yield 

(5-18) l|T(u)-r(i;)|| c([ o,7 ll Li(R*)) <CC (R)T?\\u-v\\ c[0 , nx) , 

(5.19) ||T(u) - r(«)|| C ([o,ri,£-(R*)) < CC {R)T^-i^\\u - v\\ c%n ,X)- 

Combining ([535]) . (J53JJ), (|5TT5|) and ([535]), we thus get ([5~TT|) for a G (0, 1], with 7 = | - |i, now with a new 
constant C\{R) = C(a,d,q* , RG'(2R)), where wc have chosen q* > d to ensure < |. 
As far as the continuity of T(u) with respect to time is concerned, it is enough to study 

5(u)= / Ve 5( *- s)A .*(w(s))ds. 
Jo 

We fix t € [0, T] and write for /i small enough (and positive if t — 0, negative if t = T), 

S(u(t + h)) - S(u(t)) = / Vc Ht+h~s)A . $( u ( s )) ds + / Vc 5 ('- S ) A • (e 6hA *(u(s)) - *(«(«))) ds. 

Jo 

As above, use two key estimates (|5.16l) and f|5 . together with (|5.9[) (and the dominated convergence theorem) 
to conclude the proof of Proposition 15 . 21 for a G (0, 1]. 

The case a G (1,2). We argue as before, using (|5.13[) and (|5.14[) . We need now to state and to prove the 
corresponding key technical lemma. 

Lemma 5.6. For a G (1, 2), p > and u G H a ~ 1 'P(R d ), 

(5.20) ||G(«) - G(v)\\ Ha -i, P < C 2 (\\u\\ Ha -i, P + \\v\\ Ha - 1 .,)\\u - v\\ Ha -i, P , 
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where C2(\\u\\ffa-i, P + \\v\\fja-i, P ) depends on \\u\\ H a-i, P , \\v\\ H a-i, P and on the W 1 ' 00 -norm of G' in the interval 

[Ojulloo + Hloo]. 

Proof. We use the classical identity 

(5.21) G{u) -G(v) = K{u,v){u- v) with K(u,v) = [ G'(tu + (1 - t)v) dr, 

Jo 

and the Moser estimate for the product of two functions in H a ~ 1 ' p (M. d ), see, e.g., [3SJ Ch. 2, ineq. (0.22)], to 
obtain the following inequality 

(5.22) ||G(u) - G(v)\\ Ha -i, P < \\u - «||oo||if(«,«)||B-i.p + II" ~ w||tf-i„||Jr(u > t0|| oo . 
Moreover, we recall [351 Ch. 2, Prop. 4.1] that for every increasing locally Lipschitz function H we have 

H^CiOH*.-!,, < c|if / |(|| u || 00 )|| u ||^-i,p 

for every p G (1, oo) and a — 1 G (0, 1). Choosing H = G', we deduce that 

\\K(u,v)\\ H a-l, r < f \\G'(TU+(l-T)v)\\ Ha -l, P dT 

Jo 

(5.23) < c|G"|(IM|oo + ||«||oo)(||«||b-— x., + 

Moreover, we have the trivial estimate 

(5.24) proton*, < |G'|(Nloo + Hoo). 

Combining ([5T2"2"|) . ([533"]) . ([5T2"4")) and ([33]) . we finally complete the proof of inequality ([5~2U[) . □ 

Now, we are in a position to obtain the estimate of the first component of the norm of X. From (|5.14p . we 
get 

(5.25) (u) - < C 2 {R)R\\u - v\\ p , + C 2 (2R)\\v\\ pl \\u - v\\ x < C 3 (R)\\u - v\\ x 

with C 3 (R) = CR(C 2 (R) + C 2 (2R)). We used inequality (|5T2"0|) twice, as well as the fact that ||u|| r < C||u||x 
for all r G [l,oo] and u G X. From (|5.25[) and (|5.9p with (3 = 1 and (p, r) = (1, 1), we get inequality (|5.18[) 
where Co(R) is replaced with Cq(R) = C 3 (R). 

The estimate of the 7J a_1 ' p -norm is obtained analogously. First, we have 

(5.26) - *(v)\\ p < \\u - «||ao||V a - 1 G(«)|| p + ||»||oo||V a - 1 (G(u) - G(v))\\ p < C 3 (R)\\u - 
From inequalities (|5.26[) and (|5.9p with /3 = 1 and j3 = a and (p, r) = (p, p) , we get 

(5.27) ||7» - r(u)|| c([0 ,T] > H«-^ (r-)) < CC 3 (i?)(T^ + r*-3)||u - v\\ c[0 ,t],x)- 

Finally, combining (|5 . 18[) and (|5.27[) , we complete the proof of (|5.11l) with 7 = 5 — § and with some C\ (R) . 
The time continuity of T(u) is proved as in the case a G (0, 1], and this achieves the proof of Proposition I5T21 

□ 
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Regularity of the solutions. 

Corollary 5.7 (Regularity of the solutions). Consider u G L 1 (E d ) n L°°(R d ) if a E (0,1] (as before), and 
uq G L 1 (M. d ) n (n p>Pa H a ~ 1,p (M. d )) if a £ (1,2) (a strengthened assumption). Then the solution constructed in 
Proposition \5.3\ enjoys the following regularity 

(5.28) u G C 1 ((0, T),L p (R d )) n C ((0, T), i? 1,p (M d )) 
/or every p G (p Q , oo) mi/i 

p = JjRT=J if a €(0,1], 
° l^T */<*G(l,2). 
In particular, u = u(t,x) is a weak solution of equation (|5.1|) . i.e. 

<9 t u = o"Au + V • *(n) with *(u) = |u|V a_1 G(u) 

m (0,T) x M. d in the sense of distributions (cf. equation (|5.5|l ). and 

(5.29) y u(t, x)da:= / uo(a;)da; 
for allte (0,T). 

Proof. If a G (0,1], for every u G C([0,T],Jf), where the space X is defined in (|5.10[) . we obviously have 
u G L co ((0,T),L co (R d )). We derive from inequality flU) that for all p G (p Q ,oo), we have V a_1 G(u) G 
L°°((0,T),iP(E d )) withpe = d _ ( d _ a) > 1. This implies 

V • *(«) G ((0, T), ff- 1 'P(R d )) 

for all q G (1, oo) and p G (p Q ,oo). Thus, the maximal regularity of mild solutions for the nonhomogencous 
heat equation [T7] gives us Vu G L q ((0,T), L p (R d )) for every q G (l,oo) and p € (p Q ,oo). Consequently, 
V • (MV^Gtw)) = A + f 2 with 

(5.30) /i = V|«| ■V a - 1 G( !1 ), / 2 = -|u|(-A)*G(«). 

First, we remark that /i G L 9 ((0, T), L p (R d )) for every g G (l,oo) and p G (p a ,oo). Second, we notice that 
G(u) G L q ((0,T),H^ p (R d )), hence (-A)f G(u) G L q ((0, T), H 1 ~ a > p (R d )) C L 9 ((0, T), L p (M d )). Hence, we also 
have /2 G i 9 ((0, T), L p (M. d )) for every g G (l,oo) and p G (p a ,oo). Using again the maximal regularity result, 
we obtain that 

dtueL q ((0,T),L p (R d )) 

for every q G (1, oo) and p G (p a , oo). Thus, using the following representation in L p (M. d ) for all < s < t 

u(t) — u(s) = / d t u(r) dr, 



and the Holder inequality we obtain for every p G (p a , oo) 



||u(t)-«( fl )|| p < f \\d t u(r)\\ p dr < ( f \\d t u(rW p dr) (t - s)* , 
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i.e. u £ C 0J3 ((0,T), L p (R d )) for all p £ (p a ,oo) and (3 £ (0,1). Using inequality (|5.15|) . this estimate implies 
that \&(ti) £ C°'P((0, T), L p (M. d )) for every p £ (p a ,oo). Now, the classical theory of linear parabolic equations, 
see, e.g., HH Ch. 4, Theorem 3.5], implies that d t u £ C 0,,3 ((0, T),L P (W 1 )) and u £ C°^((0, T), H 1 *(R d )) which 
is the desired regularity result. 

If a £ (1,2), we have u £ L q ((0, T), H a ~ 1 ' p (M. d )) for all q £ (l,oo) and p £ (p a ,oo), where p Q = p a = 
d/(a-l). Equations (|3.5p and (|5.20p imply that 

(5.31) V"- 1 ^) G i 9 ((0, T), L p (R d )) 

for all p £ (p a ,oo) and q £ (l,oo). In particular, *(u) = I^V" -1 ^) £ L 9 ((0, T), L p (R d )). Hence, the 
maximal regularity gives us 

(5.32) u£ L' 1 ((0,T),H 1 - p {R d )) 

for all p £ (p a , oo) and g G (1, oo). We now write once again V • = fa + fa with 

fa = sgnu Vu • V a_1 G(u), / 2 = |ti|V 0_1 (G / (u)Vu). 
In view of (|5.31|) and (|5.32p . we have 

h £ L q ((0,T),L p {R d )). 

We now claim that 

fa £ L q ((0,T), H- 1+( - 2 -^' p {R d )). 

Indeed, 

G'(u)Vu G L q ((0,T),L p {R d )). 

This implies that 

V Q_1 (G'(ii)Vu) G L q ((0,T) 7 H- 1+ ^~ a ^ p (R d )) 

which, in turn, implies the claim. Hence, 

V • *(«) G L q {{0, T), iJ- 1+ ( 2 - Q )'P(]R d )) 

for all p £ (p a , oo) and q £ (1, oo). Then, the maximal regularity implies that 

u £ L q ((0,T),H 1+ ( 2 - a) > p {R d )) 

for all for all p £ (p a , oo) and q £ (1, oo), thus we see that the space regularity of u is improved. More generally, 
the same argument shows that if 

u £ L 9 ((0, T), H l3 ' p (M. d )), with $ < a, 

then 

u £ L q ((0,T),H^ +( - 2 - a ^ p (R d )). 
Now choose the least integer k > 1 such that = 1 + k(2 — a) > a, and notice that < 2. Then 

/a G L q ((0,T),H^- a - p (R d )) c L«((0, T), L p (M d )). 
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Then the maximal regularity implies that dtu £ L q ((0, T), L p (M. d )), which implies, as was in the case a £ (0, 1], 
that u £ C°'^((0, T), L p (R d )) for all p £ (l,oo) and (3 £ (0, 1). The previous reasoning in spaces of the form 
L q ((0,T),Y) extends readily to spaces of the form C°'° ((0,T),Y). This yields the desired regularity result in 
the case a £ (1, 2). 

Moreover, it is known by [20l Ch. 4, Theorem 3.2] that mild solutions of the equation dtu = 5Au + V • / with 
/ £ C°'P((0, T), L p (R d )) are in fact weak solutions, i.e. they satisfy the equation in the sense of distributions. 
Under these regularity properties, the proof of the mass conservation property (|5.29[) is completely standard. □ 

Convexity inequalities. First, we show a simple but useful technical result involving monotone functions and 
the fractional Laplacian. 

Lemma 5.8. Let a £ (0,2]. Assume that g, h £ C^O, oo ) are strictly increasing functions. Then, for every 
nonnegative v £ C^°(R d ) we have 

h(v)(-A)^g(v)dx > 0. 



Proof. Notice that for a = 2 this lemma is obviously true, which one checks integrating by parts. Now, let 
a £ (0,2). Since (— A)^C = for every constant C £ R, we can assume that g(Q) = 0. In the same way, we 
can assume that h(0) = 0, because /(— A)%wdx = for every w £ C^°(R d ). Defining w = g(v), it suffices to 
show that 

J hig^iw^i-A^wdx > 
for all w £ C^°(R d ) such that w > 0. To do it, notice that / £ C 2 [0, oo) defined via the relation f(s) = 
/o 8 % _1 ( r )) dr for s > is convex (it suffices to check that f"(s) > 0). Hence, using the pointwise inequality 

(5.33) (-A)^g(v) < g'(v)(-A)%v, 

(see, e.g., 0, [H Lemma 1]) we obtain / h(g- 1 (w))(~A)%wdx > f(-A)%f(w)dx = 0. □ 

Next, we formulate a crucial technical tool used in the derivation of various integral estimates for solutions 
of the regularized problem (|5.ip . 

Proposition 5.9 (Convexity inequalities). Consider a C 2 function (p : R — > R + such that, for all r £ R, r ^ 0, 

<p"(r) > 0, and 

(5.34) tp(r) + \tp'(r)\ + <p"(r) < C(\r\ Ah + \r\ Ah ) 

for some constant C > and M\,Mi £ [l,oo). Then for all < s < t < T , the function u given by 
Proposition \5.S\ satisfies 

(5.35) J (p(u(t,x))dx+J J ip{u{T,x))(-A)%G(u(T,x))dxdT+S J J <p"(u)\\7u\ 2 dxdr < J (p(u(s,x)) dx, 
where tp(r) = \r\ip'(r) — f(r) sgnr. 



The proof of Proposition l579l is more or less classical, and we recall it in Appendix for the sake of completeness. 
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Remark 5.10. Remark that ip'(r) = \r\(p"(r) > 0, hence the function ip is increasing. Since G is also increasing, 
the result stated in Lemma 15.81 can be applied to show that the first dissipation term 

J iP(u(t, x))(-A)%G(u(T,x))dx dr 

is nonncgativc. The fact that this quantity is finite is a part of the result stated in Proposition 15.91 Moreover, 
PH Theorem 2.2] implies that g v (G{u)) G L 2 ((0, T), _ff f ,2 (R d )) for a function g v constructed from <p, see [IS] 
for the detailed presentation. The special case tp(r) = \r\ p is treated below. 

Remark 5.11. The convexity of <p also implies that the second dissipative term in (|5.35j) is nonnegative. Hence, 
Proposition l5.9l implies that J ip(u(t,x)) dx decreases along the flow of the regularized equation (|5.1[) . 

Corollary 5.12 (Estimates of the L p -norms). For all p G (1, oo) and < s < t, 

(5.36) J \u(t)\ p dx + {p - 1) J \u\ p - 1 u(-&)$G(u)dxdT + 8p(p- 1) J J \u\ p ~ 2 \Vu\ 2 dr < J \u(s)\ p dx. 
In particular, for p > p a (see Corollary \5. 7\ ), 

(5.37) — J \u\ p dx < -(p- 1) /" |u| p - 1 w(-A)tG(u)dx-5p(p-l) / |u| p - 2 |Vu| 2 . 
Thus, for all p G [l,oo], i/ie norm \\u(t)\\ p decreases as t increases. 

Proof. If p > 3, we can apply Proposition 15.91 with < s < t < T, </?(r) = \r\ p ; indeed, in this case, <p 
is a C 2 -function and satisfies the growth assumption (|5.34j) with {M\,M2} = (p,p — 2). Next, since u G 
C 1 ((0,T), L p (R d )), we obtain (|5.37|) from the inequality in Proposition 15.91 by a direct computation. We leave 
the details to the reader. 

If p € (1,2), we consider, for each 77 > 0, the function <p, ; such that ^(0) = f' v (0) = and 

rt(r)=p(p-l)((r 2 +r, 2 )i- l -r, p - 2 ). 

In particular, the function ip^ satisfies the assumptions of Proposition 15.91 hence, we have 

(p n (u(t,x)) dx + [ f ip r) (u(s,x))(—A)'2G(u(s,x))dxds 

<W S ,.x))|Vu( S ,,)| 2 d.xd S </^K S ,,))d.,, 

with ip' tl {r) = \r\ip'^(r), ip v (0) — 0. Letting now 77 — > and using the Fatou lemma yields the integral formulation 
of inequality (|5.37l) . 

Thus, we just proved that ||u(f)|| p < ||wo||p for all t > and p G (1, 00). By computing the limits as p — > 00 
and p — > 1, the bounds ||u(i)||i < ||wo||i and ||u(t)||oo < ||uo||oo are also obtained. □ 



We can now complete the proof of Theorem 15.11 

Proof of Theorem \5.1i In view of Proposition 15.21 and Corollary 15. 7\ it remains to prove that solutions are 
nonncgativc if initial data are so, and that solutions are global in time. 
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The positivity property is derived immediately in a usual way from the conservation of mass property (|5.6 
and the monotonicity of the L 1 -norm. Indeed, (|5.6p yields 



u(T,x)dx = j u + (T,x)dx — j u-(T,x)dx = j ((ua)+ — (uo)-)(x) dx 
and 

J \u(T,x)\dx = J u+(T,x) dx + J U-(T,x) dx < J ((uo}+ + (uo)-)(x) dx 
(here, as usual, u + — max{0, u} and u_ = max{0, — it}). These inequalities imply J U-(T, x) dx < J(uo)-(x) dx, 
and, in particular, the assumption (ito)- = gives us u > a.e. 

As far as the global existence of solutions is concerned, we argue as follows. 

For a G (0, 1], the time interval, where a solution is constructed via the Banach fixed point theorem, depends 
only on ||uo||i + || u o||oc, and this norm of the solution does not increase. Hence, we extend u = u(t,-) to the 
whole half-line [0,oo), step-by-step. 

For a G (1, 2), the Duhamel formula (|5.8|l and inequality (|5.20p with v = yield 

|Mi)||ffa-i,j> < \\u a \\ Ha -i, P + Ci(\\u \\ p *) (t- s)-%\\u(s)\\ Ha -i, P ds. 

Jo 

Due to the singular Gronwall lemma, see, e.g., [20\ Ch. 5, Lemma 6.7], we deduce that the norm ||u(t)||ij<»-i,p 
cannot explode in finite time. This shows that local-in-time solutions of the regularized equation (|5.ip can be 
also continued to global-in-time ones. □ 

6. HYPERCONTRACTIVITY AND COMPACTNESS ESTIMATES 

Hypercontractivity estimates. We now turn to prove certain L 1 H > LP estimates for solutions of problem 

no. 



Theorem 6.1 (L p -decay of solutions to the regularized problem). Let u = u(t, x) be a solution to the regularized 
problem (|5.ip constructed in Theorem \5.1[ There exists a constant C = C(d,a,m) > such that for all e > 0, 
5 > and p G [1, oo], 

d(m-l)/p+ a / V 

(6.i) \Ht)\\ p < cWuoW^- 1 ^ iwl 1 -;) 

for all t > 0. 

Proof. We first remark that it is enough to prove the decay estimate (|6.ip for large p J s, since the general result 
follows by the interpolation of the L p -norms combined with the estimate ||u(f)||i < ||uo||i from Corollary [5321 
This is the reason why we will now prove ()6.ip for p > max{m — 1, l,p a } = Pm (sec Corollarv l5 . 71 for a definition 
of p a ). 

We also remark that we can assume that M = \\uo\\i = 1 by rescaling the solution u in the following 
way. First, we consider the function u(t) = jju ( M m-i ) which satisfies equation (|5.ip with suitably rescaled 
parameters: 5 = M L-i and (if applicable) e = jr. Scaling back, we recover the desired inequality (|6.ip . 
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We first prove the result when Corollary 15.121 holds true with G(r) = \r\ m ~ 2 r and in the differential 
form (|5.37[) . This is the case when m = 2 or m > 3. In the case m £ (1,3), Corollary 15.121 holds true 
only in the integral sense (|5.36p and for a regularized function G. We will see below how to pass to the limit as 
the regularization parameter e goes to and get (|5.36[) with G(r) = r| m ~ 2 r. For expository reasons, we prefer 
to present the proof when we indeed have a differential inequality, and then to explain how to adapt it if only 
an integral version of it is available. 

The proof on Theorem 16.11 in the cases m > 3 and m = 2 is split into two steps: first, we show inequalities 
(16.11) with non-optimal constants C which blow up for p = oo; then, we improve those constants by an iteration 
method. 

Decay estimates with optimal exponents and nonoptimal constants. Our computation consists in getting the 
following differential inequality for p £ (p m , oo). 

Lemma 6.2. There exists a constant K — K(p,m) > independent of e > and 5 > such that K and K^ 1 
are bounded as p — > oo and 

(6-2) ^NI?<--K1I«II2, 



with a defined in (|3.9I) . 

Proof. We get from Corollarv l5.12l 



(6.3) J M p dx < -(p-1) J u\u\ p -\~A)%(u\u\ m - 2 )dx 



4p(p-l)(m- 1) 



a (. p + m-1 
V2 U — 2— 

(p + m - l) 2 v 1 

after applying the Stroock-Varopoulos inequality (|3.4j) with w = u\u\ m ~ 2 and q — — ^ + 1. We use next the 
Gagliardo-Nirenberg inequality from Lemma 13.21 combined with ||u(f)||i < ||uo||i = 1 to estimate the right-hand 
side of the above inequality. Thus, we get the differential inequality (|6.2|) with the constant 

4(m - l)p(p - 1) 



(6.4) K = 



C N (p + m- l) 2 ' 

□ 



With the differential inequality (|6.2p in hand, a direct computation shows that every nonnegativc solution of 
the inequality -^f(t) < —Kf(t)p has to satisfy the algebraic decay 

/(*)<(* 

We recognize (|6.1[) for p > p m with the constant C p = [K — 1^ P . Here, let us notice that the constants 
C p = C(d,a,m,p) which are obtained at this stage of the proof blow up as p — > oo. Thus, we cannot get the 
i°°-bound directly in this way. 
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Recurrence step. To improve the constant C p and to handle the limit case p = oo, we apply a variation on the 
Moser-Alikakos method of estimating the L p -norms with p = 2™ recursively, see, e.g., Q] and [TH1 Lemma 3.1]. 
The starting point is the already obtained estimate for p = 2 k > m — 1 with the least integer k. 

Lemma 6.3. For each n > k, the following estimate 

(6.5) IK*)||2» < Knt-^ for all * > 0, 

holds true with fx n = }^ n _ 1 and with a positive K n satisfying the recursive estimate 



2-(t+ i f)^+ i /(^: 



2" > 



(6.6) k„+i < 

where K n is given by formula (|6.4I) with p = 2™ +1 . 

Thus, having this estimate we see that limsupn^.^ K n < oo (irrespective of the value of at the beginning 
of the recurrence), essentially since J2^Lk n ^~ n < 00 j see Appendix IBI for details. Recall that the constants 
K = K n in the preliminary estimates have been such that K n and K~ were bounded uniformly when p = 
2™ -> oo. 

Proof of Lemma [KM We combine (|6.3j) with the Nash inequality (|3.10j) and two Holder inequalities with 2 n+1 < 
r = 2 n+1 + m - 1 < 2™+ 2 



ll 2 „ + i < IMI7IMI^ 7 with - = 2 

7 



2 n+l 



and 



As the result, we have 



ora+2 

\u\\,<\\u\\ s 2n+ 4u\\^ s with 5 = 2 . 



wwtt < -^i| U |i 2n+1 ( i+ t+^ i ),u.„- 2 "(^+^ i ) 



( j^ii-il2" + 1 - JV II "■M 2 "+i ii~iiz" 

with some if as in (|6.2[) corresponding to j> = 2" +1 . Next we estimate the L 2 "-norm of u(t) using (|6.5I) . and 
arrive to the differential inequality for f(t) = \\u(t)\\2„+i of the form 

Finally, we integrate this inequality on [0,t] and take a suitable negative power to arrive to (|6.6[) . □ 

Interpolation between p = 2" and p = 2 ,l+1 and the passage to the limit p — > oo finish the proof of the 
hypcrcontractivity estimates in the cases m > 3 and m = 2. 

To deal with the case m £ (1,3), instead of the differential inequality from Corollarv l5.12[ we use its integral 
counterpart (|5.36p . In particular, inequalities (|6.2p and (|6.7|) have their integral counterparts as well (see 
inequality (|7.5j) . below). Hence, the previous proof works in this case by applying the following lemma with 
g(t) = t and g(t) = Ct v for some well chosen positive constants C, successively. 
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Lemma 6.4. Consider functions f : [0, T] — > (0,oo) nonincreasing, and g : [0, T] — > (0, oo) increasing, smooth, 
and g(0) = 0. Assume that for a.e. t £ (0, T), s < t, 

f{t) + K J' f(rr +1 g '(r)dT<f( S ). 

Then for a.e. t G [0,T], we have 

/(*) < (JC7fl(*))"*- 

The proof of this lemma is given in Appendix C for the readers' convenience. Now, the proof of Theorem 16. II 
is complete. □ 

Compactness estimates. Now, we prove estimates which will allow us to pass to the limit as e — > and 
S — > 0, successively, in the regularized problem (|5.ip . We are going to use the fact that the approximating 
functions G = G E axe ^-Holder continuous with *y — ltiiii^tti — 1)1} on the interval [0, ||wo||oo 

I , uniformly in 

e G [0, 1]. These estimates are formulated in the following technical lemmas. 

Lemma 6.5. Assume that G is ^-Holder continuous with 7 = min{m — 1,1}. For every a G (0,1] and 
m > 1 + G (1, 2), we have 

(6.8) WV^Giu)^ < CWuW^ 

with 

q > (m - 1 - (1 - ctj/dy 1 and p = (1/q + (1 - aj/d)' 1 . 

Proof. Here, it suffices to combine the estimate on the Riesz potential with ||G(u)|[ p < C||w||^ p . The choice of 
q (or equivalently, the restriction on m) ensures that jp > 1. □ 

The compactness of a sequence of solutions to the regularized problem (|5.1[) in the case a G (1,2) is a 
consequence of the following estimate. 

Lemma 6.6. Consider a ^-Holder continuous function G : R — > R. Then, for every a G (1,2) such that 
(2 — 7)0; < 2, we have 

IIv^g^ii^cimi^ 

with 

a — 1 d a d 



7 p 2 2 

Remark 6.7. Note that the assumptions of Lemma 16.61 ensure that p > 2. 



Proof of Lemma \6.6[ We use successively the fact that H s ' p (W i ) = Fp 2 (R d ) [22j p. 14], the characterization of 
Triebel-Lizorkin spaces with difference quotients 1221 p. 41], and finally known embedding theorems for Besov 
and Triebel-Lizorkin spaces [221 p. 31] in order to derive 

llV^GHHp < CHG^OH^-i < C[GUu\\ ^± < C\\u\\ H§>2 , 

P,2 p 7 11 
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where [G] 7 := sup r .^. s \G(r) — G(s)\/\r — s| 7 and p > 2 is chosen so that the so-called differential dimension is 
constant; this yields the condition appearing in the statement of the lemma. The proof is now complete. □ 

Lemma 6.8. For a £ (1, 2) and (3 — m)a < 2, there exists p £ (1, oo) such that 

(6.9) ||V a - 1 G(«)||,<C||«|| jr i. a(Rd) 

where C depends on the "/-Holder seminorm of G in (0, ||u||oo) {with 7 = min{m — 1, 1}). 

Proof. Apply Lemma 16.61 with 7 = m — 1. □ 

Lemma 6.9. If a £ (1, 2) and m > a, there exists p £ (1, 00) and r £ (1, 00) {and r > m — 1) such that 

(6.10) ||V a - 1 |«r- 1 s g a«|| P <C|||«| £t ¥ =1 sgii«|| ff9 , a(R(l) . 

Proof. Consider v = \u\^^ sgnu. Then estimate (|6.10[) is equivalent to the following one 

\\V^\v\^sgnv\\ p <C\\v\\ Hf ^ Rdy 

We then apply Lemma |6~61 with 7 = ^^^rj arL d get the desired result if 7 > 2(1 — ^), or equivalently, 

m — 1 

r < . 

a- I 

Here, in order to find r £ (l, ""j 1 ), the condition m > a is needed. The proof is now complete. □ 

7. Proof of Theorem 12.61 

This section is devoted to the proof of our main result on the existence of solutions to problem Ql.ip — (|1.2p 
satisfying decay estimates ([2.60 . 



Proof of Theorem \2.b\ We first consider very regular initial data, i.e. we assume that uq satisfies (|5 ,3[) . Then, 
this condition is relaxed by considering initial data that are merely intcgrablc. 

The proof proceeds in three steps: passage to the limit with the parameter of the regularization of the 
nonlincarity, then with the parameter of the parabolic regularization, and finally — stability with respect to 
initial data. 

Passage to the limit as e — > 0. Consider «o satisfying (|5.3p . From Theorem 15.11 we have a sequence of 
solutions u £ of (|5.ip for G e defined in such a way that there exists C > such that for all e £ (0, 1], 

[GA :=SU P ' G f -y>l < 
r^s \r-sp 

where 7 = min{m — 1,1}. Thanks to Corollary 15. 12[ there exists also a constant C > (depending on S > 0) 
such that for all e £ (0, 1], 

(7-1) H u elll,~((UT)xR £! ) < G, ||'" e ||j,2((o,T),.H" I . 2 (R<i)) < G. 

Hence, we can construct a sequence {e n )„ such that 

U n -± u in ^((O.T),^ 1 '^^)) 
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where u n denotes u £n . Moreover, for all R > 0, the embedding H 1 ' 2 (Br) c L 2 (Br) is dense and compact. 
Using (|7.ip we also obtain 



lim / / \u n (t,x)\ 2 dtdx = 0. 

meas(-E)->0,.EC[0,T] J £ J Bh 

Hence, we infer from |21j (which contains an optimal result on the compactness for Hilbert space valued vector 
functions) that, up to a subsequence, for every R > 0, 

u n ^um L 2 ((0,T) x B R ). 

Moreover, passing again to a subsequence if necessary, we can assume that for every R > 0, 

(7.2) Vu n Vu in L 2 ((0,T) x Br), u n ->■ it for a.e. (i,a;) G Q T , 

which imply for all 93 £ C^°((5t), 

/ / u n dt^p dt dx — > / / w9tv?dida;, 
JJot JJqt 



Vu n ■ dt da; — > / / Vu • V<^ dt da;. 

J Jq t 

As far as the nonlinear term in equation (|5.1j) is concerned, we have 

|w„ — > \u\ for a.e. (t,x) € Qt- 
Hence, the dominated convergence theorem combined with (|7.1|) imply that 

\u n \V(p^ \u\Vip in L 2 ((0,T),L q '(R d )) 
for all g' 6 (l f oo). Lemmas 16.51 and [6T8l imply that 

(7.3) V a_1 G £n (u n ) is bounded in L q (Q T ) 

for some g 6 (1, 00). 

We deduce from (|7.3p that we can extract a weakly converging subsequence of (\7 a ~ 1 G 6n (u„))„ C L 9 (Qt) 
which limit is denoted by L. Since 

G En (u n ) — > |it| m_1 sgnu a.e. in Qt, 

we have 

V a ~ 1 G en (u n ) — L in L q (Q T ), where X = V a_1 (|u| TO - 1 sgnu) in V(Q T ). 
In particular, V Q_1 (|u| m_1 sgnw) £ L q (Qr) holds, and 

V a - x G en (u n )-\u n \V(pdtdx^ [[ V^dur^sgnuJ-lulV^dtda; 

J Jq t 

as n — > 00. Hence, we obtain a weak solution of (|5.ip with G(r) = rj"^ 1 sgnr for uq satisfying (|5.3|) . 
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By the Fatou lemma, we derive from Corollary 15.121 and Proposition 13.11 the following estimates: for a.e. 
t > s > 0, q S [1, oo], p € (1, oo), uniformly in 5 € (0, 1], 

(7-4) K*)ll,<l|uo||„ 

(7.5) N*)llg+ 4P r (P 7 1)(m i : 2 1) // vHsgnulul^) 2 dsd X <\\u(sW p , 

(7.6) ||«(t)||, < C\ | u 1 1 f^S? r ^t^W (i - i ) , 

(7.7) 2(5 // |Vw| 2 dsda; < ||u 



112- 

Remark that now we can derive the hypcrcontractivity estimates in the case m £ (1, 3) from inequality (|7.5p . 
We also know that the solution w is nonnegative if uq is so. 

Passage to the limit as 6 — > 0. Let u = u s denote the solution constructed above. We consider 

v s = \u\ m sgnu. 

Using (|7.5p with s = and p = m + 1 > 1, we get that, for all T > 0, 

v s is bounded in L 2 ((0,T), H^' 2 (R d )). 
Hence, there exists a subsequence S„ — > such that 

« n -w in.L 2 ((0,T),iff' 2 (M d )), 
where v n denotes v dn . Moreover, inequality (|7.4[) with q = oo implies that for i? > 

lim / / \v n (t,x)\ 2 dtdx = 0. 

mcas(£;)^0,£C[0,T] ^ 7 Br 

Hence, we can use [5T| once more, and conclude that for all R > 

w n ^v in L 2 ((0,T) x 
Moreover, passing to a subsequence if necessary, we can assume that 

u n — > u for a.e. (t, x) € Qt, 
where u n denotes u Sn . Now Lemmas 16.51 and 16.91 imply that 
(7.8) V^^Kp" 1 sgnuj is bounded in L q (Q T ) 

for some q > 1. Hence, we can pass to the limit in the nonlinear term in the weak formulation of equation (|5.ip . 
To complete this proof, notice that inequality (|7.7[) implies that 

S n Vu n -> in L 2 (Q T ). 

In particular, 



<S n / / Vu„ • V<p di dx — s- 
as n — >• oo. We thus conclude that u is a weak solution of (|l.lj) - (|1.2|) . 
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The conservation of mass, the positivity property, the monotonicity of L p -norms, and the hypercontractivity 
estimates follow from (|7.4[) - (|7.7|) in a standard way. 



Stability with respect initial conditions. Assume now that u g i 1 (R d ) fl L°°(M. d ). Consider an approx- 
imating sequence Uq satisfying (|5 . 3[) . Then, the sequences of i p -norms are bounded and we can pass to the 
limit as we did already when letting <5 — ^ 0. We also recover the expected properties of the weak solution. 
Assume finally that u is merely integrable. Then T n (u ) G L 1 (M d ) n L°°(R d ) where 

T n {i~) = min{max{r, — n}, n}. 

Hence, there exists a weak solution u n of (jl.ip . In view of the hypercontractivity estimates, we can extract 
a converging subsequence in (rj,T) x M. d , with arbitrary r\ > 0, to a function u v in the following sense 

u n v? in L 2 ((7],T) x B(0,R)). 

Moreover, the L 1 -norm of u n is bounded in (0,T) x M. d . Then we have 



/ / u n dt<pdtdx 

ii(0,i))xl d 



<C V , 



u n d t (p ->• / / u^dt^ 

(r;,T)xR d J J {■q,T)xM d 



sgnu„) • V<^di da; -> // V c 

V7(, ; ,T)xR d 

as n — > oo for each r\ > 0. We can now conclude the proof of Theorem 12.61 through a diagonal procedure. □ 



V a ~ 1 (|u n r -1 sgiiu„)-Vpdfda;-> // V a_1 (|u 1j r -1 sgnu") • V^dt da; 

(r;,T)xR d 



Appendix A. Proof of Proposition 15.91 
Proof of Provosition \5.9\ We fix an arbitrary p > p a = — ^j. First, we remark that 

J J {d t u)(j)dtdx + J HV a-1 G(u) ■ V<t>dtdx + S J J Vti- \J(j)&t&x = 

for all compactly supported in time and in L 1 ((0,T),H 1 'V' (R d )). Recall that we have d t u G C[0, T], i p (M d )) 
for p > p Q = ^ and = |u| V Q_1 G(u) £ L°°((0, T), L p (M. d )). To justify the previous equality, it is enough 
to mollify the function (f> in time and space, and remark that the mollified function ^ satisfies (for r\ small 
enough) <f> n {0,x) = <f> v (T,x) = 0, u<f> n G C 1 ((0, T), L 1 (K d )) and 

u(dt<f>n) dt dx = — J J (dtu)(f> v dt dx. 

Letting r\ — > and using the regularity of ty(u) yields the desired result. 

Next, consider <fi(r,x) — ip'(u(r, a;))0„(r), where is truncation function in time of [s,tj. Remark that 
4> e L 1 ((0, T),H 1 ' P ' (R d )); indeed, (p'{u) £ L°°((0, T), L p '(M d )), and we also have for all iu £ (H 1 ^nL 1 nL 00 )(R d ) 

\\<p"(w)Vw\\pi < \\Vw\\ p \\if"{w)\\ q 
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with i = 1 — -. As far as 9,, is concerned, we choose it such that Q v (s) = and O^(t) = p^ij— s— rj) — p^ir— t+rj) 
where p v is an even mollifier supported in [—rj, rj\. Now we can write 



(d t u)ip' (u)Q v drdx = jj <9 t (<y9(it))8 r) dr da; 

^(w)(/o I? (r - t + rj) — pr,(T - s - rj)) drdx. 



Hence, 



Moreover, 



{dtujtp' {u)Q v dr dx — > / f(u(t, x)) dx — / <f(u(s, x)) dx. 



|u|V<0(u) = \u\<p"(u)Vu = V(tp(u)) in L p (R d ) 
thanks to the Stampacchia theorem (u <S L°°(M. d ) hence ip' and ip locally Lipschitz is enough). Hence 

j J |tt|V a-1 G(u)-V^(tt)dTda: = J J V a_1 G(u) • V(tp{u)) dxdr. 
It remains to prove that 

(A.l) / ( iP(u)(-k)^(G(u))dTdx < [ /v a_1 <3(«)-V(V>(u))da;d<r<oo. 



The last inequality comes from the computations we made above. As far as the second inequality is concerned, 
we use (11.31) to write 

dx dy dr 



J J V^Giu) ■ V(V(u)) dx dr = Jim j* j j {G{u){t, x) - G(u)(r, y))F v (y - x) ■ V(^(u(r, x))) 



where F v is defined as follows 



F v (z) = C d 



Through an integration by parts, we now get 

(G(u)(t, x) - G(u)(t, y))F n (y - x) ■ V{^(u(t, x))) dx dy dr 

VG(u)(t,x) ■ F n (y — x)4'(u(t,x)) dxdydr 

+ / (u(t,x) -u(T,y))(V • F v )(y -x)ip(u(T,x))dxdydi 



The first term on the right-hand side equals since F v is odd. Moreover, we have 

r 1 d +« + (d + a + l)\z\ d + a 

Hence, the Fatou lemma yields (|A.1[) . The proof of the proposition is now complete. □ 

Appendix B. From (|6.6[) to the boundedness of k„ 
Consider l n = log n n and write (|6.6p as follows 
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with 



On+l (aim 
Z V d + 2" I 

0„ = 1 



TTTT log 



2™ 


' 2a 




+ l\ 




id + 2" 




; 



Remark next that 



In particular 



f 2 n + 1 + 2(m - 1) ' 



n C 

a n < C — and o„ < 1 . 

™ —2 n 2" 



a„ < oo and J 6 n < oo 



Using the fact that b n < 1, we get 



°fco- 



Hence, l n does not blow up, and neither does n n . 



Appendix C. Proof of Lemma \6A 



Proof of Lemma \6.4\ First, we remark that we can reduce to the case g(t) = t and K = 1 through a change of 
variables. 

The proof is simple if / is smooth. If / is not, extend / by to R and consider a mollifier p e . Then write 
for ti <t 2 , 

f(h -s)+K ^ P+\t) dr < f(t 2 - s). 
Now integrate against p e (s) and nse the Jensen inequality to get 

fe(h)+K r n +1 {r)Ar<f e {t2). 

J — oo 

We are now reduced to the case f = f e , which is smooth. Passing to the limit, the proof is now complete. □ 
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